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The Astrodynamical Space Test of Relativity using Optical Devices (ASTROD) is a 
multi-purpose relativity mission concept. ASTROD's scientific goals are the measurement 
of relativistic and solar system parameters to unprecedented precision, and the detection 
and observation of low-frequency gravitational waves to frequencies down to 5 X 10"'' 
Hz. To accomplish its goals, ASTROD will employ a constellation of drag-free satellites, 
aiming for a residual acceleration noise of (0.3-l)x 10-15 g-2 Hz-i/2 at 0.1 mHz. 
Noise sources and strategies for improving present acceleration noise levels are reported. 
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1. Introduction 

The classical concept of a drag-free satellit^ consists of a small proof mass inserted 
inside a large spacecraft. The larger spacecraft shields external forces, allowing 
the proof mass to move in free fall. The relative position and orientation of the 
proof mass with respect to the main spacecraft is sensed along its trajectory. This 
information is feedback to thrusters on the main spacecraft, which are subsequently 
fired to maintain the proof mass-spacecraft relative position and orientation. In this 
way, the coupling of external forces to the proof mass is minimized. 

Numerous space missions have already employed drag-free technology. The first 
example of a "drag- free" mission to test fundamental physics is Gravity Probe B 
(gp-b).I21gp-b was launched in April 2004 and has experimentally measured frame- 
dragging and geodetic effects. By spring of 2007, the results of the data analysis will 
become public. The GP-B inertial sensor achieved a level of free fall below 2 x 10"^^ 
m s-2 at 5 X lO^^ Hz. 

The Laser Interferometer Space Antenna, LISA, is another fundamental physics 
mission requiring drag-free performance. The LISA mission conceplP^ consists of 
three spacecraft in heliocentric orbits, forming a nearly equilateral triangle forma- 
tion of side 5 x 10^ km and inclined with respect to the ecliptic by 60°. LISA will 
monitor the separation between free falling proof masses, which are shielded within 
the spacecraft, by using interferometric techniques, to detect and observe gravita- 
tional waves. LISA'S aims include studying the role of massive black holes in galaxy 
evolution, testing relativistic gravity, determining the population of ultra-compact 
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galactic binaries, probing the physics of the early universe, observing supermas- 
sive and intermediate black holes mergers and mapping spacetime by observing 
gravitational captures. The LISA drag-free performance goal is 3 x 10^^^ m s~^ 
at 0.1 mHz. By 2009, the LISA Technology Package (LTP) on board the 
LISA Pathfinder (LPF) ESA mission, with NASA contributions, aims to demon- 
strate drag-free performance to a level one order of magnitude lower than the LISA 
requirement, approximately 3 x 10^^^ m s^^ Hz^^/^ in the frequency bandwidth 
between 1 mHz and 30 mHz.Hllf LTP is successful, it will achieve the best drag-free 
performance up to present date. 

The Astrodynarnical Space Test of Relativity using Optical Devices, ASTROD, 
is a mission conceplP^ that consists of a constellation of drag-free spacecraft em- 
ploying laser interferometric techniques with Earth orbiting satellites, to provide 
high precision measurements of the relativistic parameters {7, /?}; improved deter- 
mination of the orbits of major asteroids; measurement of solar angular momentum 
via the Lense-Thirring effect and the detection of low-frequency gravitational waves 
and solar oscillations. ASTROD aims to improve on the LISA drag-fre e goal by a fac- 
tor of between 3 and 10, i.e. (0.3-1) x lO^i^ ni s^^ Hz-^/^ 0.1 mHz.l3SEIllIlini 
It is worth noting that a ten-fold acceleration noise improvement with respect to 
LISA, would allow ASTROD to explore relativistic gravity to an uncertainty level of 
1 ppb F I 8 | 13 | 14 | A simple version of ASTROD, ASTROD I, has been studied as the 
first step to ASTROD. ASTROD I concept consists of one spacecraft in a solar orbit, 
carrying out interferometric ranging and pulse ranging with ground stations .^^^^^21 
ASTROD I also aims to measure relativistic and solar system parameters and test 
fundamental laws of spacetime although with less precision than ASTROD. 

Other mission concepts to follow on from LISA are (Big Bang Ob- 

server) and DECIGCM^ (DECihertz Interferometry Gravitational wave Observa- 
tory). These missions aim to observe the cosmic gravitational wave background 
produced by standard inflation, having optimum sensitivity around 0.1 Hz, where 
white dwarf binaries confusion level is thought to be very low. The drag-free force 
noise requirement for BBO and DECIGO is approximately a hundredth of the LISA 
force noise target. 

In section 2 we report general ideas that ASTROD could adopt to improve drag- 
free performance, we summarize acceleration and sensor back action disturbances, 
emphasizing the most significant low-frequency noise sources. Charging disturbances 
and discharging schemes are briefly discussed in section 3. Finally in section 4 we 
discuss the performance and problems associated with inertial sensors for other 
follow on LISA missions. Gravitational interactions are described in terms of mul- 
tipole moments and a detailed calculation of capacitances for the case of a three 
dimensional capacitive sensor/actuator are given in appendix A and B, respectively. 
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2. ASTROD Inertial Sensor 

ASTROD will face new challenges in drag-free technology compared with LISA. 
First, ASTROD aims to improve the LISA drag-free performance target by a factor 
of between 3 and 10 at 0.1 mHz. On the other hand, ASTROD will extend the 
gravitational wave observational bandwidth to frequencies below 0.1 mHz, which is 
the lowest frequency on the LISA observational bandwidth, down to 5 x 10^^ Hz. 
To achieve these goals, there are key problems that need to be evaluated and their 
solutions optimized. 

The initial design concept proposed for LISA was that laser beams from re- 
mote spacecraft would directly illuminate the proof masses. Issues like laser beam 
pointing, actuation forces to correctly orientate the proof mass with respect to laser 
beam, and cross coupling of proof mass degrees of freedom, make this scheme less 
atractive. A new scheme has been recently discussed, in which the laser beam from 
the remote spacecraft illuminates a fiducial point in the inertial sensor or spacecraft, 
rather than the proof mass. The inertial sensor provides positioning and attitude ref- 
erence, and is also used to actuate the proof mass along the other degrees of freedom. 
The proof mass could be monitored by employing heterodyne laser metrology. T his 
scheme, so-called separate interferometry, has been adopted recently by LISA.ini 
However LISA will not fully exploit all its advantages as it may still use two proof 
masses per spacecraft as references for two different interferometer arms. ASTROD 
will employ one proof mass, minimizing disturbances and simplifying control. 

Replacing capacitive sensing with optical sensing, for drag-free missions follow- 
ing LISA, has also been widely debated. Optical sensing is more sensitive than 
electrostatic sensing and it requires almost no coupling between the proof mass and 
surroundings. Towards that directio n many efforts to implement optical sensing 
have been made in the last few years ,^1^2^^21^1122] further laboratory research 
is needed to develop a space qualified optical sensing scheme. Light pressure could 
also be used for active control. A more conservative design would combine an opti- 
cal sensor and capacitive active control.^^Sl Larger gaps between the proof mass and 
electrodes could then be used to minimize disturbances. 

Ultimately local gravitational gradients are the limiting factor for drag-free per- 
formance. In that context the influence of proof mass geometries (spherical, cylin- 
drical, cubic, polyhedral, etc) on sensitivity and its repercussion for the overall in- 
ertial sensor design merits further discussion. For monitoring and correcting length 
changes due to thermal effects and slow relaxations, ASTROD will use an absolute 
metrology system.f^D 

2.1. Acceleration noise sources 

To discuss the acceleration noise, we consider a simplified control lo op model of a 
spacecraft and a single proof mass. The acceleration noise is given by,l25l26l27] 




(1) 
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Table 1. Direct acceleration disturbances. The parameters are 
defined as follows: A and E denote cosmic ray impact rate and 
incident energy, P and T are housing pressure and temperature, 
respectively, Ap the proof mass cross section, and are the 
electrostatic and magnetic shielding factors, respectively, STqb 
describes optical bench temperature fluctuations and ^yg the 
thermal shielding factor between optical bench and the proof 
mass, STsc represents spacecraft temperature fluctuations and 
a thermal expansion coefficient. Finally kg, a and G denote 
the Boltzman, Stefan-Boltzman and Newton Gravitational con- 
stants. The rest of parameters are defined in the text. 

Environmental disturbances 



Cosmic rays 

Residual gas 
Magnetic susceptibility I 
Magnetic susceptibility II 
Permanent magnetic moment 
Lorentz I 
Lorentz II 
Radiometer effect 
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Thermal radiation pressure 
Gravity Gradients 
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where Fstr are stray forces directly acting on the proof mass of mass vrip^ Fns and 
TNt are external forces and thruster force noise acting on the outer spacecraft of 
mass M . Xnr is the sensor readout sensitivity, u is the control loop gain and uj = 
27r/, where / is the frequency. External forces, thruster noise and readout sensitivity 
contribute to acceleration noise because of the proof mass-spacecraft coupling K. 
If the proof mass and spacecraft are highly decoupled, then the acceleration noise 
will be given by stray forces directly acting on the proof mass. 

In Table 1 direct acceleration noise sources with the exception of sensor back 
action disturbances are listed.^^SI Environmental disturbances can be divided into 
different groups depending on their origin. There are disturbances caused by im- 
pacts. Cosmic rays which penetrate the spacecraft shielding and residual gas can 
deposit momentum onto the proof mass. There are disturbances of magnetic ori- 
gin. Proof mass magnetic susceptibility, Xy and residual permanent moment, Mr, 
can interact with the residual local and/or the interplanetary magnetic field, Bsc 
and Bip, respectively. There are also disturbances associated with charge. Residual 
charge accrued on the proof mass can interact with the interplanetary magnetic 
field via Lorentz force. Finally, there are disturbances associated with thermal fluc- 
tuations on the spacecraft. These include the radiometer and outgassing effects, 
thermal radiation pressure and gravitational gradients caused by thermally induced 
spacecraft distortions. 
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Several key factors needed to quantify the total acceleration noise, such as sensor 
readout noise, back action forces, and stiffness terms, differ for different types of 
sensor. As an example, optical sensing is nearly stiffness free, providing high readout 
sensitivity with very low back action forces. On the other hand, with the widely used 
capacitive sensor, high sensitivity is achieved at the expense of increased acceleration 
noise and stiffness. Also, because of the fact that close metallic surfaces are needed 
(3-4 mm gaps in the case of LISA), other noise contributions due to, for example, 
patch effects and dielectric losses, become significant. 

In the best scenario, the proof mass will ultimately be coupled to the space- 
craft by gravitational gradients. Force gradients of electrostatic origin can be made 
negligible by implementing large gaps between the proof mass and the surrounding 
metallic surfaces. An analysis of proof mass geometries and thermal-gravitational 
modeling of the spacecraft and payload are necessary to account for gravitational 
disturbances and stiffness terms. 

A preliminary quantitative analysis of acceleration noise parameter requirements 
for ASTROD is given in Refs. Eland [H 

2.2. Low-frequency acceleration noise sources 

The LISA observational bandwidth extends from 0.1 niHz to 0.1 Hz. It has been 
pointed out that gravitational wave observations extended to frequencies below 0.1 
niHz, are desirable in the study of certain astrophysical sources like massive black 
holes (MBH) binaries at high redshift.^^ 

The ASTROD free falling proof masses will be separated by distances of 30 to 
60 times longer than those of LISA. ASTROD gravitational wave sensitivity curve 
will therefore be shifted to lower frequencies than the target LISA bandwidth. At 
low frequencies, spurious forces acting directly on the proof mass are the dominant 
source of noise. This is the reason why, for a mission like ASTROD, it is particularly 
important to identify these sources of noise. An extended discussion of low- frequency 
acceleration noise sources and low-frequency sensitivity curve for gravitational waves 
for ASTROD is given in Ref. [H 

Thermal, magnetic and electrostatic effects are sources of low-frequency accel- 
eration noise. Thermal noise arises due to radiometer effect; fluctuating outgassing 
and thermal radiation pressure assymetries; thermal distortion of the spacecraft and 
residual gas impacts. The magnitude of these effects for a particular mission are de- 
pendent on the mission orbit, which dictates the thermal environment. Suppressing 

on 

thermal disturbances requires thermal diagnostics,*^ stable electronics, passive and, 
in some cases, active thermal isolation, thermally conductive electrodes (for the case 
of electrostatic sensing/actuation) and high vacuum. A preliminary evaluation of 
thermal disturbances for ASTROD shows that the outgassing effect, thermal radi- 
ation pressure and thermally induced gravity gradients could be at levels of about 
fog w 1.1 X 10^1^ m s-2 Hz-i/2, ftp - 8 X IQ-^^ m s^^ 11^-1/2 ^nd fgg w 5.4 x 
10"^^ m s^^ Hz^^/^, respectively, at 0.1 mHz. A vacuum pressure of the order of 
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10^® Pa, and a thermal isolation factor, ^ts, of about 150, were assumed for these 
estimates. Below 0.1 mHz, solar irradiance fluctuations become the main cause of 
temperature fluctuations. Solar irradiance fluctuations become more acute when 
approaching solar rotational period, which is of the order of 25 days.^^ This as- 
pect will be a crucial factor for the thermal diagnostic and thermal isolation system 
design for ASTROD. 

Magnetic noise at low frequencies is caused by interplanetary magnetic field 
fluctuations, local gradients and magnetic field fluctuations, eddy current damping, 
magnetic impurities, Lorentz forces due to proof mass residual charge, etc. Suppres- 
sion of magnetic disturbances requires further reduction and shielding of permanent 
magnets in the payload, improving magnetic shielding, adopting a magnetic clean 
wiring (i.e, solar array rewiring, etc) and power system. The most significant mag- 
netic low-frequency noise source is due to the interaction of proof mass magnetic 
susceptibility with the interplanetary field fluctuations, fm2 (see Table 1). Assum- 
ing parameter values given in Ref . [11] fm2 ~ 2 x 10"^^ m s~^ Hz~^/^ at 0.1 mHz. 
Disturbance fm2 increases at low frequencies as f~'^^^. 

Other relevant noise sources at low frequencies are caused by electrostatic ef- 
fects, i.e., due to voltage noise, charge fluctuations, DC voltages, dielectric losses, 
actuation noise, etc. When employing capacitive sensing, the strategies to suppress 
these noise sources are: increasing the separation between the proof mass and the 
electrodes, active compensation of DC voltages,!^ avoidance of DC voltages applied 
to drag-free degrees of freedom, high quality surface coatings to minimize dielectric 
losses, high stability power supplies and continuous discharging of the proof mass. 
Electrostatic noise is caused by the capacitive sensor /actuator. 

An obvious way for suppression of these noise sources is to replace capacitive 
sensing by optical sensing. Nevertheless, other concerns will arise if optical sensing 
is to be used. Thermal distortion of optical components and changes in refractive 
index with temperature modify optical paths. Because these noise sources are due to 
thermal fluctuations, they will also be of signiflcance at low frequencies. New ideas 
addressing these problems, such as using all reflective optics, by using gratings, has 
been extensively discussed in the literature.^1^ Diffractive gratings have also been 
considered to enhance a ngu lar sensitivity, compared with standard angular sensors 

lOOl 

based on laser reflection.'^ 
2.3. Gravitational modelling 

Drag-free performance will ultimately be limited by local gravitational fields and 
field gradients. Because of structural distortion of the spacecraft due to thermal 
fluctuations, thermal and gravitational disturbances need to be modelled together. 
For the present discussion, we will be concerned only with the local gravitational 
interaction between the test mass and a simplified spacecraft structure. 

We consider three different proof mass geometries: spherical, cylindrical and 
cuboid. For simplicity, the spacecraft will be assumed to be a hollow cylinder, as a 
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first approximation. 

Following appendix A, we can analyze the gravitational interaction by means of 
inner, qim, and outer multipole moments, Qirm that describe the proof mass, and 
the spacecraft and payload mass distribution, respectively. Using expression lA.3[ the 
first non-zero outer multipole contribution of a hollow cylinder is Q2o- Expression 
IA.5I shows that Q20 couples to the proof mass inner multipoles goO: Qi±i, <l20, 92±i 
and q2±2- 

The first non-zero inner multipole moments for a parallelepiped proof mass 
of sides 2a, 2b and 2c are goo = ^p/V^, 920 = 2/3-\/5/47rmp (2c^ — — 5^) 
and q2±2 = l/12-\/l5/27rmp(a^ — b^). On the other hand for a cylindrical proof 
mass of radius R and semi-height h, we have qoo = mp/^/47T and 520 = 
mpy/b/A^ {h? /2> — R^/A). Given these values, it can be seen that the quadrupole 
moments vanish for a cubic proof mass (as the one adopted by LISA). In the case 
of ASTROD I, a proof mass of dimensions 50 x 50 x 35 mm is considered. In that 
case q2±2 vanishes but not 920- To first order in the gravitational interaction with 
a cylindrical tube the (720 term appears as a constant energy term and will not 
contribute to the gravitational force. Disturbances proportional to proof mass cross 
section area can then be minimized by shortening one of the dimensions (as is the 
case of ASTROD I proof mass) . A trade off between the acceleration disturbances 
which are proportional to proof mass cross sectional area, and the gravitational 
interaction needs to be done. If a cylindrical proof mass is utilized we can also sup- 
press quadrupole gravitational interaction by choosing h/R — V3/2. In the same 
way disturbances proportional to cross sectional area can be suppressed by short- 
ening /i, without altering gravitational interaction with "far away" gravitational 
asymmetries. A preliminary analysis of gravitational force gradients for ASTROD 
I is given in Ref . [511 

Another issue that needs to be considered for ASTROD is the fact that the 
relative distances and angles between the spacecraft are not constant along their 
orbits. Therefore, telescopes utilized for laser beam pointing need to be steered 
during the missio n. Fo r active gravitational compensation ASTROD will employ 
dummy telescopes .1^2^ 

2.4. Back action disturbances 

When deciding which type of sensor to use, there are two main options to consider. 
First, we could consider a low-stiffness sensor. In that case, the proof mass is highly 
decoupled from the spacecraft, at the expense of losing sensitivity. The other option 
is to employ a high-stiffness sensor to achieve better readout sensitivity, at the 
expense of a high level of coupling between the proof mass and the surrounding 
structure. Capacitive sensing exemplifies this issue. To improve sensitivity we need 
to place the electrodes closer to the proof mass. By doing that, the stiffness and 
back action disturbances increase (see Table 2). 

To understand the principle behind stiffness and back action disturbances due 
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to electrostatic sensing and actuation, we first consider the total mechanical energy 
for a capacitive sensor-proof mass system. Following Refs. [25] and [27l the total 
mechanical energy is given by 

i 

where q is the net charge of the proof mass; C is the coefficient of capacitance of the 
proof mass: C = '^^Ci, where i = xi, X2,yi,y2, zi, Z2, g defines the capacitances 
formed by an electrode facing a proof mass side, and capacitance to ground, Cg; 
Vs is the voltage induced on the proof mass due to the applied voltages, Vi, and 
voltage to ground, Vg-. Vs = J2i ^i^i- 

The force acting along a generic direction, assuming neither charge and voltage 
gradients, is given by 

F=\Y.C'^{V^-VMf (3) 

i 

where Vm = + c' ^^"^ ^'i derivative of capacitance Ci along the generic 

direction. 

Force disturbances can then be considered of two types: a) position dependent 
disturbances, caused by fluctuating position and attitude of the proof mass (stiffness 
terms), and b) position independent disturbances, caused by voltage and charge 
fluctuations. Along a generic sensitive (drag-free) axis, we can write the fluctuating 
force terms due to charge and voltage fluctuations as 

5Fsv.. 1 = ( V q (y, - ) - c I sv. (4) 



5Fsv,2 ^\\^C[:^ + yU\^ a (5) 

- i i / 

5Fs,,,^V:^^5q (6) 

5Fs,.2 = V:^Jq (7) 



where 



i 

and 

i 

Table 2 shows the back action disturbances given above for the special case 
of a one dimensional capacitive sensor and one translational degree of freedom.^^SI 
Disturbances due to dielectric losses and patch flelds are also listed. An extended 
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Table 2. Sensor back action acceleration disturbances. The pa- 
rameters are defined as follows: 5i;die! denotes voltage fluctuation 
due to dielectric losses, Vo is dc bias voltage, Vpe average patch 
potential, VxO average potential across opposite side of sensor, SVj^ 
fluctuations in voltage different across opposite side of sensor and 
d is the gap between the proof mass and surrounding housing and 
electrodes. The rest of parameters are defined in the text. 

Back action disturbances 

Dielectric losses foL = VoSvaiei 

Patch fields (Uncompensated) fpe = ^ ^ ^fj-VpeSq 



discussion of electrostatic back action disturbances for ASTROD I and ASTROD 
can be found in Ref . [571 and [Til respectively. 

Position dependent disturbances (stiffness terms) can be obtained by calculating 
the variations in capacitances, SCi, and capacitance gradients, SC'^. 

In appendix B formulae for capacitances, capacitance gradients and their fluc- 
tuations, given by Ci, C-, 6Ci and (5C-, respectively, are obtained for the special 
case of a 6-degree of freedom, cubic, capacitive sensor. These expressions codify 
cross-coupling effects between translational and rotational degrees of freedom. 

On the other hand, optical sensing offers advantages in terms of high sensitivity 
and low back action forces. Optical readout sensitivity is ultimately limited by laser 
shot noise. Laser shot noise is proportional to where P is the laser power. 

Lasing back action force is proportional to the laser power, and is given by 2P/c, 
where c is the speed of light. This force can, in principle, be compensated to a high 
degree of accuracy. 

3. Charging disturbances 

Galactic cosmic rays (GCR) and solar energetic particles (SEP) incident on the 
spacecraft will result in the accumulation of charge on the proof mass. Charge 
accrued on the proof mass leads to numerous sources of noise. Some of these noise 
sources have been described and discussed above. First, when summarizing direct 
acceleration noise sources, the Lorentz force due to the movement of the charged 
proof mass through the interplanetary magnetic field was discussed. We have also 
discussed, in the particular case of employing capacitive sensing/actuation, how the 
proof mass charge couples to sensing/actuation voltages to induce spurious forces 
and stiffness terms that will affect the performance of the inertial sensor. 

Disturbances associated with charging can be divided into three types: a) those 
which are proportional to charge accrued, q, b) those proportional to Sq, which are 
so-called "shot noise" terms and c) mixed terms proportional to qSq. This division 



SVd X q 
Sq X Vd 
5q X q 



fsv.i = ^^{V,o-V,)5V, 
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is of importance in understanding disturbance suppression schemes (see discussion 
below) . 

Time dependent forces also contribute to the spectral noise density. Proof mass 
charging is a time dependent process and both Coulomb and Lorentz forces give 
rise to coherent Fourier signals (CHS). Assuming a linear increase of proof mass 
charge with time, the total charge can be written as q{t) = qt + Sq, where q is the 
mean charging rate. Following Ref. [32|the acceleration noise terms due to Coulomb 
and Lorentz interactions are given by 

acHS = hkit) ^ (0fc + Qk)t + Ekt^ (10) 

where 

0fc = T^, 0fe = TTl^d, and = — Ad (11) 

TOp^e nipC a nip \Ca J 

using a parallel plate approximation to estimate capacitances and capacitance 
derivatives. The parameters used above are defined as follows: v is the orbital veloc- 
ity of the proof mass, Bip is the interplanetary magnetic field, is the electrostatic 
shielding factor, Cx is the capacitance along the sensitive axis, Vd is the voltage 
difference between opposite sensor sides, and d and Ad are the capacitance gap and 
gap asymmetry, respectively. 

Inspection of pip shows that geometrical and electrostatic asymmetries in the 
inertial sensor contribute to the appearance of these signals. Geometrical asymme- 
try arises due to limited machining accuracy and it is represented by an asymmetry 
in the capacitance gap. Electrostatic asymmetry is due to a stray DC potential 
imbalance between opposite sides of the sensor. These residual DC stray potentials 
are dependent on the work function of the metallic surfaces. These potentials are 
measurable in average, for each electrode, and can be balanced by appropriate ap- 
plied bias volt ages .'^Ultimately, voltage offset compensation will depend on voltage 
measurement precision required and work function domains stability in periods of 
time comparable with the measurement integration time. In the context of LISA 
mission, these coherent Four ier signals can exceed the instrumental noise target, for 
typical parameter values.^^ 

At low frequencies, charging disturbances and coherent charging signals are of 
particular concern. Those so-called "shot noise" charging disturbances and coher- 
ent signals scale, roughly speaking, with frequency as 1//.-32- Of special note is the 
acceleration disturbance proportional to residual voltage difference across opposite 
sides of the sensor and charge fluctuations, fsq.i (see Table 2). By active compen- 
sation, the potential difference across opposite sensor sides can be balanced to '--^l 
mV.^f^ To improve on the LISA acceleration noise by a factor of 10 at 0.1 mHz, 

1 9 

ASTROD would need to reduce this value to below 0.7 mV, giving a disturbance 
level fsg^i w 1.4 x lO^^^ m s^^ Hz^i/^^ 
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3.1. Discharging schemes 

When discussing discharging schemes we have to keep in mind that not only charge 
accr ued by the proof mass but also the charging rate are potential causes of 
j-^Qjgg 133134] Discharging the proof mass will suppress some charging disturbances 
and will reduce the coupling of the proof mass with its surroundings. Coherent sig- 
nals are proportional to mean charging rate. By "mean charging rate" we mean that 
charging and discharging rates are added linearly to give a net rate. Therefore by ac- 
curately matching charging and discharging rates, these signals can be suppressed.^ 
Nevertheless the fact that proof mass discharging involves the transport of charge 
"packages", will cause additional shot noise. Shot noise due to proof mass charg- 
ing and discharging cannot partially cancel each other because the charging and 
discharging processes are statistically independent. Therefore their respective shot 
noise terms have to be added quadratically. 

For LISA, a continuous discharging scheme will be adopted. The continuous 
discharging process will consist basically of two steps. Firstly, the proof mass charge 
needs to be accurately measured. That can be done by applying a sinusoidal dither 
voltage and measuring the displacement along a non drag-free axis. This proof mass 
displacement is proportional to proof mass charge. Secondly, UV light will shine on 
the proof mass an d/or surro unding electrodes to discharge the proof mass via the 
photoelectric effect .1221321341 

The ASTROD strategy to suppress charging noise will depend on which sens- 
ing/actuation device is employed. ASTROD could benefit from the replacement of 
capacitive sensing by optical sensing. Even if a capacitive scheme is employed for 
force actuation, charging disturbances and coherent signals could be suppressed by 
increasing the gaps between the test mass and surrounding surfaces. If optical force 
actuation is employed, then only Lorentz type disturbances need to be considered. 
In that case, the charging requirements can be relaxed and the discharging scheme 
can be simplified. 



4. Inertial sensors for other missions to test fundamental physics 

Follow on LISA mission concepts have been proposed, not only to extend the ob- 
servational bandwidth to lower frequencies, but also to fill the gap between space 
antennae and ground based interferometric facilities (LIGO, GEO600, VIRGO, Ad- 
vanced LIGO, LCGT, etc). Follow-on mission concepts include the Big Bang Ob- 
server (BBO) and the Japanese antenna DECi-hertz Interferometer Gravitational 
wave Observer (DECIGO). Both BBO and DECIGO are to be designed to have an 
optimum sensitivity between 0.1 Hz and 10 Hz. The objectives of the decihertz an- 
tennae are: measuring the expansion rate of the universe; determining the equation 
of state of dark energy, by observing the coalescence of binary neutron stars and 
stellar mass black holes; and shedding light on the growth of superm assive black 
holes, by studying the merger of intermediate mass black holes! ^ ^ 

These missions will also measure the relative distance, or maintain the distance 
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by feedback, between test masses in nearly free fall. To achieve their objectives they 
need to reduce residual spurious forces to approximately a hundredth of the LISA 
goal. 

These missions are conceptually different. To shift the gravitational wave sen- 
sitivity curve towards the decihertz level, the effective interferometer arm length 
has to be, approximately, a hundredth of the LISA arm length (5 x 10^ km). The 
BBO preliminary conceptual design consists of a constellation of spacecraft with 
a LISA-type design and an arm length of 5 x 10^ km. On the other hand, DE- 
CIGO plans to place in space a Fabry-Perot cavity of length 1000 km and finesse 
of approximately 10. To accomplish this, DECIGO would place massive mirrors of 
100 kg mass and 1 meter diameter in free fall. Large actuation forces, to keep the 
cavity in resonance, will be required. This condition and the stringent residual ac- 
celeration noise levels, of the order of 4 x lO^^'' m s^-^ Hz^^/^, seem to be difficult 
to reconcile. Nevertheless one could employ a large control loop gain to minimize 
the acceleration disturbances due to actuation forces. Other technological difhcul- 
ties common to both missions are related to the requirement for low residual gas 
pressure. Extremely low pressure could be achieved by venting some of the gas to 
outer space. However, this could cause other problems such as drag of the proof 
mass because of residual gas flow, and undesired particles coming from the thruster 
propellant, brought into the proof mass housing. 
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Appendix A. Gravitational interaction in terms of multipole 
moments 

The gravitational potential can be written in terms of multipole moments aJ^^^Sl 

oo / ^ 

V^-AttGY, E ^f-Y<llmQlm (A.l) 
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where qim (Qim) are the inner (outer) moments defined respectively by 

qim= [ pt{x'y%lAe',cj,')<fx' (A.2) 

J Vt 

and 

Qhn^ f p.(f)r-('+')ri„(0,0)d3f (A.3) 

J Vs 

If the inner multipoles of the test mass, qim, are known in a given reference frame, 
then the inner multipoles with respect to a new reference frame, can be obtained. 
Let us consider a position r" — r' + f. The test mass position with respect to the 
reference frame, of origin O, in which the multipoles are known are denoted by f. 
On the other hand, r* denotes the position vector of the origin O, with respect to 
the new reference frame in which we wish to work out the multipoles. In the special 
case of pure translations,!^^! 



;j'=Om,m' V ' J \ 'I 

5l,i+i'C{1', m\ I, m, L, M)Yi1„^,{e' , mrn{e, </)) (A.4) 

Using eq. IA.4I we can rewrite the gravitational potential energy lA.ll in terms 
of the known inner multipoles and the translational parameters between the two 
reference frames as 



LI' —0 m,m' 



The gravitational force can be then easily calculated by taking the derivatives 
of lA. 51 with respect to the translational parameters. 



Appendix B. A cubical inertial sensor. Capacitance calculations 

We consider a six-dimensional degree of freedom model for the capacitive sens- 
ing/actuation device. A cuboid proof mass of side lengths {2Lx, 2Ly, 2Lz) is inserted 
into a three dimensional capacitive sensor. The gaps at the equilibrium position be- 
tween the proof mass and the electrodes are denoted by {D^, Dy, D^). The proof 
mass translational degrees of freedom are denoted by {dx,dy,dz). The proof mass 
rotational degrees of freedom are given by the Euler angles (0, 0, in the so-called 
" x-convention" . 

We first consider the electrostatic energy density to work out capacitances for 
this model. If the electric field between two conducting surfaces is defined by i?, then 
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the electrostatic energy density is given by a; = ^EqE^. The proof mass faces and the 
electrodes will define the capacitances. In the case in which the displacement and 
rotation of the proof mass are infinitesimal, we can approximate the electric field 
between conducting surfaces by Ei ~ where Vm (Vi) and Axi denote 

the proof mass (electrode) potential and the capacitance gap in the i-direction, 
respectively. The electrostatic energy is then given by 



W^{d,,dy,d,,^,9,i^) ^ -eoiVM -V,f / — ^ (B.l) 



where by the simbol ± we differentiate between the gaps at opposite sides of the 
sensor. 

By integrating along the gap (in this case we choose the gap along the z-axis), 
the electrostatic energy can be written as 

W^id,,dy,d^,cb,e,^)^-leoiVM-V,r I '!'"^]> , „ (B.2) 

2 J Ax- (x,y,d^,dy,d^,(j),0,%p) 

and we can define the capacitances by 

Ct,{d,x,dy,d^,(j),9,ip) ~ eo [ dxjdxk—^ ^ i n i\ (^-3) 

J Ax^ [Xj,Xk,dj:,dy,dz,(p,0,^p) 

The information about how different degrees of freedom couple to each other is 
codified in IB. 31 

To explicitly work out capacitances, we translate and rotate the proof mass by 
the parameters (d^, dy,dz, (j), 0, ijj). 

The rotation matrix in terms of the Euler angles (0, 0, ip) is written as 

cos tp cos (j) — cos 9 sin </> sin if] cos sin + cos cos </> sin ip sin sin \ 
- sin -0 cos 4> — cos 9 sin (f> cos ip — sin ip sin (j) + cos 9 cos (p cos 4' cos ip sin 9 (B.4) 
sin 9 sin (j) — sin 9 cos (j) cos 9 j 

Given this matrix we can define the normal vectors to the proof mass faces after 
a ((/), 0, '0)-rotation, with respect to a fixed reference frame. This reference frame has 
its origin in the geometrical center of the sensor and its axis orthogonal to electrode 
surfaces. The normal vectors can be written as 

' cos ip cos 4) — cos 9 sin </> sin ip ^ 
cos V' sin + cos cos (/) sin | (B.5) 
sin sin 9 

■ sin cos — cos 9 sin <p> cos '0 ^ 
sin sin + cos cos cos -0 | (B.6) 
cos ipi sin 9 

























2) \ 




1^3 = risy ~ — sin(?cos0 (B.7) 
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The equations of the proof mass faces, are defined by 

(x-P^)-n, = (B.8) 

where Pi = ±Lini + d, being d = (d^, dy, dz) a displacement, and Li defines the 
semi- length of the proof mass along the three axes x, y and z. 

Using this expressions we can obtain the capacitance gaps for the three direc- 
tions. These are given by 

Ax^ ^F.^±{y + At)n^,±{z + Bt)n^2 (B.9) 

Ay^ - F± ±{x + A^) n\ ± (z + i3±) n| (B.IO) 

Az^ = Ff ±{x^ A^)n\±{y + B^)nl (B.ll) 
We define the capacitances as 

Gtu^ =eo [ ' dy [ ' dz_J__ (B.12) 







-L, Aa;±(?;,z) 



/-L^ 



Ly 



1 



Cfdo^n - £0 / dy dz—rr—. (B.13) 



Ax^{y,z) 



Ct.,H. = / . i^l dy ■ (B.17) 



pL^ ^ 



-Ly Az^{x,y) 



where by up, down, left, right, we indicate that two electrodes face each proof mass 
face (see Fig. [T]). 

The capacitances in the x-direction are given by 

C+ - ^° £ In + + ^^^'"^ + + ^"^"^ 
"^"f 712 " F+ + {B+ - Lz)n^^ + {A+ +Ly)nl 

eo F+ + (L, + B+)n- + Atnf F+ + (L, + B+)n- + {A+ + Ly)nl 



F+ + (L, + B+)n^ + A+nf 



so F+ + {Bt - Lz)n- + Atnf F+ + {Bj - L,)^! + {A+ + L,)nf 



ASTROD Inertial Sensor 17 



c 



'y.down,- 



Fig. 1. Schematic view of capacitive sensing electrodes. 



^{L, + B-)nl-{A- +Ly)nl 



In- 



F." - {B^ - L,)nl - {A-^ + Ly)nl 
F- -{L, + B-)nl-{A- +Ly)nl 



Fx - {Lz + Bx )nf - Ax nf 
-{B- -Lz)n-^~{A-+Ly)nl 



Fx - {Bx - L,)nl - Ax 



(B.19) 



+ ^£0. , F+ + + Bt)ni + (A+ - 

" F+ + {Bt -L,)n- + {A+ -Ly)nl 



n. 



so F+ + (L, + B+)n- + Atn- F+ + (L, + + (A+ - 



F^ + [L, + B+)nl + A+nf 



eo F+ + {B+ - L,)n^ + f + + - L,)n| + (A+ - L,)nf 



^- ^_£ir 1 - (-^^ + Bx )nl - {A~ ~ Ly)nf 

.,dorun ^x y p- _ _ L^)n- - (Ax ^ Ly)n- 

en F- -{L, + B-)nl-A-nl F' - {L, + B-)nl - [A- - Ly)nl 

nf nf Fx - (L;: + Bx)n^ - Axnf 

gp F- - {B- ^ L,)nf - A-nf - {B- - L,)n^ - {A- - Ly)nf 

nf nf Fx - {Bx - L,)n^- Axnf 



where 
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F^ = L,{l-n,,) + D,Td, (B.22) 

At = TL^niy - dy (B.23) 

Bt = TL.ni, - 4 (B.24) 

< = ^ (B.25) 
nix 

n| = — (B.26) 
nix 



The capacitances in the y-direction are 



^y,up — ^x,upii^v^ ^z, A^ ,111,712) {Lz,Lx,By ,Ay ,n|,n^))(B.27) 



Ctao^n = Ctio^^{{Ly, L,,AtBt,nl, nl) (L„ Lx, B^, A±, n|, n^))(B.28) 



where 



=Ly{l-n2y) + DyTdy (B.29) 
A^ = TLyU^x - dx (B.30) 
B± = TLyn2, - d, (B.31) 

(B.32) 



n2x 

n2y 

y _ 



nl=— (B.33) 
ri2y 



The capacitances in the z-direction are given by 



^f,ieft ~ '^t,up{Lz ^ Ly-, Ly Lx) (B.34) 

^f, right = ^x,down{^z Ly', Ly Lx) (B.35) 



where now 



Ft = L,{1 - 713,) +D,Td, (B.36) 

Af = TL,n3x - dx (B.37) 

Bf = TL.nsy - dy (B.38) 

= ^ (B.39) 

n'2 = — (B.40) 

We can approximate the expressions of capacitances by taking into account that 
F > A and B. Then we have 

i=7 + -Bjnf + [AZ + Ly)nf 
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To work out the force and the force disturbance along the sensitive axis we also 
need to work out capacitance gradients, which are given by 

(B.45) 
(B.46) 
(B.47) 
(B.48) 



r)r<+ 


2eaLyLz 


ddx 


[F+ + B+n^ + {At + Ly)nl] ^ 




2eoLyLz 


ddx 


[Fx -Bxnl^{Ax +Ly)nlY 


x,down 


2£t^LyLz 


ddx 


[F+ + B+nf + {At - Ly)nlf 


dC~ 

Xydown 


2£oLyLz 


ddx 


[Fx - BxTi^ - {Ax - Ly)nfY 


Along the y-axis we have 




2eaLxL,X4> + "0) 


ddx 


[F+ +A+n\ + {B+ +L,)nlY 


dC~ 
y,up 


2eoLxLz{(f) + 


ddx 


[Fy- -Ay n\- {By +1^)711]" 


dC+^ 

y,aown 


2eoLxLA(P + ^) 


ddx 


[Fy++A+n\ + {Bt-L,)nl]' 


y^down 


2eoLxL,{cf) + tjj) 


ddx 


[Fy--Ayn\-{By-L,)nl]' 


Along the z-axis, 






2eoLxLy(j)0 


ddx 


[F+ + B+n'2 + iAt + Lx)nl]^ 




2eoLxLy(j)6 


ddx 


[F^ - B^n'2 - i^z + Lx)nl]^ 




2eoLxLy(p9 


ddx 


' [F+ + B+n^^ + {At -Lx)n(\^ 



(B.49) 
(B.50) 
(B.51) 
(B.52) 

(B.53) 
(B.54) 

(B.55) 
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^^z,right 2soL^Ly(j)0 



The terms useful for stiffness calculations are variations of capacitances and 
capacitance gradients. These are given by 

^Ct^, = . ^ 1^ ^ [^-^ + ^^^"^ + {At + Ly)n^,] (B.57) 

[F+ + B+nl + {At + Ly)n\] 

5C-u^ = 6 [F- - B-nl - {A- + Ly)nf] (B.58) 

[Fx - B:r nf - [Ax - Ly)nf\ 

Finally the last useful formulae for stiffness calculations are those of the type 
5C-. On the x-axis we have 



[F+ + B+nl + {At + Ly)nl] 

- \F,.Bt^ti-.y,r^ ' - ^^^"^ - ^''^ - '^^^ 

<.o^n = ^.^"ff' ,,..3 ^ [^x- - B-n% - {A- - 

[Fx -Bxul- {Ax - Ly)nl\ 

On the y-axis the expressions are as follows: 
^'"^ [Fy+ + Atnl + {Bt + L.)nlf ^ ' 

[F+ + A+nf + {B+ - L^)nl\ 

xn-' ^SoLxLz{4> + 1p) r y T \ VI 



And finally on the z-axis we have 



dC. 



z,left 



P^^^^^4^---^S [F/ + + {At + Lx)nl] 



[F+ + B+nl + {At + Lx)nl 



B.61) 
B.62) 
B.63) 
B.64) 

B.65) 
B.66) 
B.67) 

B.68) 

B.69) 
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where for the x, y and z axes we have 

S [F+ + Btnl + {At ± Ly)nl\ = -6d, ± L,(^<^ + 5i>) (B.73) 

(5 [F- - B-n^ - (A- ± Ly)nl\ = +5d^ T Ly{6<p + 5il)) (B.74) 

(5 [F+ + A+n\ + {B+ ± L,)n^] = -6dy ± (B.75) 

d [F- - A-n\ - [B- ± L,)nl] = ddy T LJ0 (B.76) 

d [F+ + Bfn^ + (At ± L^)n(\ = -Sd^ (B.77) 

S [F- - B-n'^ - {A- ± L,)nf] = 5d, (B.78) 

Appendix B.l. Capacitance as a position sensor. 

By measuring and combining capacitances along the different axis, we can obtain 
the position and attitude of the proof mass. The parameters {dx,dy, dz, 4>, 0, ip) are 
obtained by the following combination of capacitances, 

- C-J + (C+ - C-,) SeoL^Ly^ (B.79) 

(C+r - C-r) - (C+ - C-i) Seo%^0<^ (B.80) 



{(^x,up - Cx,up) + i^x,down ~ '^x^own) - ^^O^zLyj^ (B.81) 
iC+,up - C-,^p) - iC+do^on - CxAo^n) ^ -^SoLIL, + i^)^ (B.82) 



T, T. 

' Dl 



{^y,up - ^y,up) + i^y,down ~ '^y,down) — ^^oLxLz^ (B.83) 



iCy,up - Cy^up) - {Cy^down " ^y,down) - -^SoL^L10^ (B.84) 
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Appendix B.2. Capacitances, capacitance derivatives and their 
variations 

For the special case in which the proof mass is in the equihbrium position, 0, with 
no translational and rotational offsets, 



ri± — r<^ _ 

^V,up ^y,down p, 
z,left z,right 

ddx ddx 

9Cy,up _ ^^'y.down _ ^^'z'jeft _ ^^'f,right ^ q 

ddx ddx ddx ddx 

^C^,up = T^%^ [-6d. + Ly{5<t> + 6n 

^x 

= \-Mx - L,{64> + Sn 

SC^,up = T^^^[-Sdy + L.5e] 

V 

<do^n = T^-^^[-Sdy-LJ9] 

y 

5Ct^p = ^^^^ [-5d^ + Ly{6<f> + 6^,)] 



5C: 



x,down ^3 



^^y,up — ^^y,down — ^^z,left " ^^z,right ~ ^ 



